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We consider fluctuations of the electronic spin due to coupling to nuclear spin. Noise spectroscopy 
of an electronic spin can be revealed in the Scanning Tunnelling Microscope (STM). We argue that 
the noise spectroscopy of electronic spin can reveal the nuclear spin dynamics due to hyperfine 
coupling. Tunnelling current develops satellites of the main lines at Larmor frequency and at zero 
frequency due to hyperfine coupling. We also address the role of the rf field that is at or near the 
resonance with the nuclear hyperfine field. This approach is similar to Electron Nuclear Double 
Resonance (ENDOR), in that is allows one to detect nuclear spin dynamics indirectly through its 
effect on electronic spin. 
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I. INTRODUCTION 

Noise spectroscopy is a technique that allows one to 
measure spectroscopic properties by observing nontrivial 
features in the noise. The main feature of noise spec- 
troscopy is that the noise spectrum has encoded in it 
spectroscopic features that correspond to physical exci- 
tations in the system, such as atomic levels, Zeeman split 
levels of electrons or nuclear levels. The transition be- 
tween levels would cause enhanced dissipation when the 
energy transferred equals the energy difference between 
those levels. This enhanced dissipation also would imply 
enhanced fluctuations at the same frequencies/energies, 
as follows from the fluctuation-dissipation theorem. Ex- 
periments that prove utility of the noise are available 
from many fields. An incomplete list includes nuclear 
spin noise, 1 Faraday rotation noise in the alkali atoms, 2-4 
and acoustic noise. 5 Recently noise spectroscopy was 
used to detect a single electronic spin with Magnetic Res- 
onance Force Microscopy by an IBM group. 6 

Therefore, in principle, noise measurement could be a 
powerful tool to investigate the dynamics of the system. 
Sometimes noise is easier to measure and then noise spec- 
troscopy could be even a preferred technique to investi- 
gate nano-scale systems. 3,6 

One example of the noise spectroscopy relevant for us 
here is Electronic Spin Resonance Scanning Tunneling 
Microscopy, ESR-STM. ESR-STM is a technique that is 
using the extremely local nature of the STM measure- 
ment to detect the noisy precession of spin centers on the 
nonmagnetic surface. When a tip of an STM is located 
above a paramagnetic spin center the tunneling current is 
modulated by the precession in the presence of external 
field. It was shown 7,8 that the ac current at the Larmor 
frequency is spatially localized within 0.5 — lnm. It is 
the spatial localization that suggests that this technique 
is capable of detecting a single spin. In addition it was 
proved that the frequency of the signal is dependent on 



real time on the size of the magnetic field. 9,10 More re- 
cently similar experiments have been done on the param- 
agnetic BDPA molecule. 11 The interest in this technique 
has risen sharply recently, due to the possibility to ma- 
nipulate and detect a single spin 12,13 and due to the pos- 
sibility to use it for quantum computation. 12,14 There 
have been many proposals for the mechanism of this 
phenomenon. 13,15-20 Present experiments are not suffi- 
cient to constraint possible mechanism and further inves- 
tigation would help to elucidate the nature of the effect. 

Recently Durkan reported the measurements of the 
noise in the ESR-STM on a TEMPO molecule. 21 This 
molecule is well characterized, and it contains nitrogen N 
with the nuclear spin 1=1. ESR spectrum of TEMPO 
in the bulk is known to exhibit the hyperfine splitting 
on the order of 15 G that corresponds to the free 
(<7 e = 2) electron precession frequency on the order of 
45 Mhz. The main new observation that is important in 
our context is that the ESR-STM on TEMPO molecule 
produced three peaks that possibly corresponds to the 
hyperfine split Larmor line in current spectrum. If re- 
produced, this observation opens up a new possibilities 
in noise spectroscopy in STM. 

The purpose of this paper is to consider the case of 
coupling of the electronic impurity spin S to nuclear spin 
I via hyperfine coupling. We investigate the noise spec- 
troscopy of the nuclear spin and coupled electron-nuclear 
dynamics as it might be seen in STM experiments. In- 
deed all the previous discussions of the mechanisms so far 
have focused on the dynamics of the electronic impurity 
spin. Hyperfine coupling to the nuclear spin would allow 
one to measure dynamics, relaxation times and transi- 
tions between the nuclear levels. 

We find that the hyperfine coupling produces addi- 
tional satellite lines in the noise spectra for the localized 
impurity spin that will be split away from Larmor line 
by the amount proportional to hyperfine coupling A. We 
also consider the case of the rf field at frequency ui ~ A 
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applied to the system. This is the case of Electronic- 
Nuclead Double Resonance (ENDOR), where electron 
spin dynamics will be affected by the nuclear spin flips 
caused by the rf field. In the case of applied rf field we 
clearly pump energy into the system and hence measure- 
ment is not strictly a noise spectroscopy measurement. 
Still this is a set up most likely to be attempted experi- 
mentally and this is why we address it here as well. 

The plan of the paper is as follows. We address the 
localized impurity spin susceptibility with the hypcrfinc 
coupling and the tunneling current modulations in the 
Sec. II. We address a few specific cases, such as case of no 
rf field and case of no applied fields, neither dc nor ac. In 
all of these cases electronic spin will have nontrivial spec- 
troscopic features, most notable the hyperfine split lines 
around Larmor line and around zero frequency. In Sec 
III we conclude with discussion of possible experiments. 
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FIG. 2: The possible spin-flip transitions that may occur in 
the system. The electronic spin-flips give rise to the three 
major peaks in the transport current fluctuations spectrum 
at the frequencies E z and E z ± 2A; see also Fig. 3 and Fig. 
4. 



II. PROBING THE SPIN SUSCEPTIBILITY VIA 
NOISE MEASUREMENT 

To be specific, consider the tunneling current between 
two contacts in the presence of a localised spin S inter- 
acting with a nuclear spin I, / = 1, to make a contact 
with the experiments on TEMPO. The Hamiltonian of 
this system is written in the form 

H = H l +Hr + H s + J2 4j*o + *iS-a a , a ]c p/ 3, (1) 

kpaft 

where H L[R) = T,k( P )a e K P )acl {p)lJ c k(p)a models free 
electrons in the left (right), L(R), lead, whereas 

Us = B a S z + AS -1 + B X I X cos cut + B 2 S X cos ut, (2) 
5 > A^>B ll B 2 

accounts for the interactions between the localised elec- 
tronic and nuclear spins. In the Hamiltonian, Eq. (I), 
a a p denote the Pauli spin matrix vector with matrix in- 
dices a, (3, the Fermionic cL \ a i c k(p)a- are creation and 




FIG. 1: Sketch of the experimental setup. The tip — surface 
separation distance is d. The single magnetic impurity atom 
with spin S on the substrate surface interacts with the nuclear 
spin I in the substrate. 



annihilation operators of electrons in the fcth (pth) eigen- 
state in the STM tip (substrate surface), with spin a =f 
,|. For simplicity we have incorporated the electronic 
g e and nuclear gyromagnetic ratio into the effective 
fields B = g e H , B\ = gNHi,B 2 = g e H\ for the exter- 
nal dc field Hq and ac field Hi cos cut. The last term in 
Eq. (1) describes tunneling of the electrons from the tip 
into the substrate surface in the presence of the localised 
electronic spin S. This term only give a contribution to 
the net steady state current by providing a chemical po- 
tential shift (bias voltage drop) between the two bands 
{£ka\ and {Spa}- In real systems, the hopping matrix 
elements t = to + iiS • a have a kp— dependence, which 
are omitted here in order to make the notation more com- 
pact. The wavefunctions of our system are superpositions 
of the direct product states I^l) <S> |V"s) ® (V'fl)j e -g- the 
direct product of the state of the STM tip, the impurity 
spin, and the substrate surface. The tunneling matrix t 
in the last term of Eq. (1) couples all of these different 
states, of which the term proportional to to describes the 
spin independent tunneling while the term proportional 
to t\ provides the spin dependent contributions arising 
from the exchange interaction of tunneling electrons to 
the magnetic atom. 

Spin dependence of the tunneling arise due to di- 
rect exchange dependence of the tunnel barrier. 22 The 
overlap of the electronic wavefunctions of the tip and 
the surface, separated by a disctance d is exponen- 
tially small and is given by a spin dependent tunnel- 
ing matrix element i = 7exp{— a/(<E> — JS(t) ■ ct)/$o}, 
where the direct exchange between the tunneling elec- 
tron spin a and the impurity spin S is explicitly in- 
cluded. Here, J is the exchange interaction parame- 
ter between the electrons tunneling from the tip to the 
surface and the precessing impurity spin S. The tun- 
neling barrier height $ is typically a few eV, whereas 
$ = h 2 /(8md 2 ) is related to the distance between the 
tip and the surface. 23 Since the exchange term in the ex- 
ponent is small compared to the barrier height, we may 
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expand the exponent in JS, which gives t = to + tiS(t) • 
a. Here, to = 7 cxp (-y/®J$o) cosh [JS/(2y/$$ )] de- 
scribes th e spin independent tunneling, while t\ = 
7exp (-a/^/^q) sinh [JS/(2y/^ )} accounts for the 
spin dependent tunneling amplitude. For estimates we 
may employ the typical thumb rule ii/to ~ JS/{2§) <C 
1. 

We will assume that the tunneling electrons are par- 
tially spin polarised. This can be achieved in several 
situations. For example, in ferromagnetically coated 
tips a potential difference 2Sfi a separates the spin bands 
£ka = £k + uSficr- Fcrromagntically ordered tips have 
proven successful in the study of magnetic structures. 24 
Another approach to obtain a spin polarised current is to 
use an antiferromagntically coated tip with no ferromag- 
netic order. 25 Such tips have the benefit of a vanishing 
dipolar field and should, therefore, have a negligible influ- 
ence on the precession frequency of the impurity spin. A 
ferromagnetically ordered tip may produce a field of 0(1) 
T at a separation of few Angstroms from the surface. 
Such fields therefore leads to huge precession frequencies 
which are difficult to measure. In what follows, however, 
we will not be concerned of how the spin polarised current 
is generated. We define a parameter A which relates the 
spin polarised current to the net tunneling current. This 
is the parameter that will be determined by a particular 
microscopic model of the tip. Thus, we will henceforth 
treat A as a phenomcnological parameter. 

For conciseness, we employ the Heisenberg picture for 
all operators {O}, i.e. 0(t) = exp [iHt]Os exp [— iHt], 
with Os = 0(t) for operators in the Schrodinger picture. 
All finite temperature expectation values (0(t)) will rep- 
resent T,iPi(M0)\O{t)\tpi{0)), where ^(t = 0) is the 
zero time wave function from the Schrodinger picture, 
whereas pi is its probability within the density matrix 
formulation. 

Now, for a qualitative description of the effect 
addressed here consider the charge current J(t) = 
—ed/dtNL{t). Since we are considering the steady state 
regime it is, by charge and current conservation, sufficient 
to consider the current in the tip (or in the substrate) 
only. By a direct calculation, using the Heisenberg equa- 
tion of motion, we find that 

J(t) = -2elm ( c Lfa> + tiS(t) • a a(3 ]c p/3 ), (3) 

kpafi 

where e is the electronic charge. Hence, we see that the 



tunneling current can be partitioned into a spin inde- 
pendent part J (t) = -2elm^2 kpcr to (c^Cp^) and a spin 
dependent part 

6J(t) = (SJ(t)) = et 1 (S(t)-3 s (t)), (4) 

which depends on the localised moment S(t), where 

J s (t) =-i^2 C ka a aP C P0 + H c -> ( 5 ) 

kpaft 

is the spin dependent contribution to the tunneling cur- 
rent. The z-componcnt of this expression, e.g. Jg{t) = 

~i Sfcp( c It c pT — c ki c pl) + H.c, describes the net flow of 
spin f and spin J. carriers, whereas the transversal com- 
ponent (Jf, jv)(t) = -tJ2k P ( c U c Pl + c li c Ph ~ lc U c Pl + 
icLcpf) + H.c. accounts for spin flip transitions of the 
tunneling electrons caused by the interactions with the 
precessing impurity spin. 

To the lowest order in the tunneling amplitude t±, 
the electronic current-current correlation function arising 
due to the spin dependent part of the current is given by 

<{<JJ(t), JJ(f)}> = (ei 1 ) 2 (5 i (i)^(t'))(JJ(i)Ji'(t')) 

+(*<->*'). (6) 

where ({•,•}) denotes the symmetrized correlator whereas 
i,j = x,y,z signify the spin components. Thus, to low- 
est non-trivial order in t\ we can treat the two tempo- 
ral correlation functions K^(t - t') = {S l {t)S 3 {t 1 )) and 
C{t-t>) = (Ji(t)Ji(t>)) - (Ji(t))(Ji(t')), (\t-t'\ ^ 00) 
independently. To make a connection with the main pro- 
posal of this paper we note that the Fourier transform of 
the symmetrized correlation function ({63(t), SJ(t')}) is 
the current noise spectrum at various frequencies arising 
from the localised electronic spin. In Fourier space, the 
current noise spectrum is given by the convolution of the 
two power spectra associated with S and a, e.g. 

( |SJ (£ )|2) = I K v(s')C(e - e')de> + (e - -s). 

(7) 

In order to see the effect of the interactions between 
the localised electronic (S) and nuclear (I) spin in the 
rotating magnetic field we have to calculate (e) and 
K z {e), for which the details are given in the appendix, 
see Eqs. (A6) and (A8). The resulting expressions for 
the correlation functions are given by 
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In the above equations we have defined the detuning 
parameter Alu = A — uj/2 and the parameter wi = 
\J Aui 2 + (Bi/2) 2 . Resonant conditions of the system is 
given for Aio = u> = 2A giving u>\ = B\/2. While 
the coefficients (B\l2uj\) 2n n = 1,2, rapidly decays to 
zero out of resonance, the coefficients (Aw/wi) 2 ™, n = 
1,2, rapidly grows to unity, since the amplitude of the 
rf field B\ <C A. Hence, all terms in Eqs. (8) and (9), 
but the ones proportional to cos(w ± 2wi)r/2, contribute 
to the spectrum at resonance, whereas only the peaks at 
e = 2Bo and e = 2Bq ± (u ± 2u>) give a non-negligible 
contribution to the spectrum out of resonance. Below we 
will focus on few experimentally relevant possibilities. 



1. B =£0, Si ^ 

This case corresponds to applying an external DC 
field and rf field in or close to nuclear resonance, 
which is known as Electronic Nuclear double Resonance 
(ENDOR). 26 Under those circumstances the nuclear spin 
dynamics is influenced by the rf field, which is reflected 
in the dynamics of the electonic spin. This is perhaps 
one of the most relevant cases for nuclear spin noise to 



be seen in STM experiments. 

Clearly, the correlation function K + ~ in Eq. (8) has a 
central peak at the Larmor frequency e = 2Bq, which at 
resonance is represented by the first, third and seventh 
terms in Eq. (8). This peak should also be measurable 
out of resonance since the constant (first) term is non- 
vanishing for all frequencies of the rf field. There are 
side-bands around the frequency e — 2Bq ± 2A, repre- 
sented by the second, third, fourth, tenth, and eleventh 
terms, of which only the fourth and fifth are present out 
of resonance. Hence, also these peaks should be seen for 
all rf frequencies. However, at resonant conditions there 
are peaks around e = 2B ± A, represented by the sixth, 
eighth, ninth, twelfth, and thirteenth terms in Eq. (8), 
which are not expected to be measurable appreciably far 
out of resonance. The same observations hold for the 
correlation function K z , although its spectrum is cen- 
tred around e — 0. Hence, the central peak of this part 
is hidden in the white noise spectrum. 

In Fig. 3 we display a contour plot of the correlation 
functions K + ~ (upper panel) and K z (lower panel) as 
function of the Fourier frequency e and detuning Aw. 
Especially for it is readily seen that there are five 

peaks (2_B , 2_B ± A, 2B ± 2 A) around resonant con- 
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FIG. 3: Contour plots of the spectral intensities (e) (up- 
per panel) and K z (e) (lower panel) as function of the detun- 
ing Auj in Fourier space (e). The spectra are generated by 
Lorentzian functions of uniform width T/Bo = 1/60, using 
A/ B = 1/12 and Bi/B = 1/100. 

ditions Aoj = 0, while only three peaks remain out of 
resonance. In Fig. 4 (solid) we show that total spectrum 
at resonance for the same conditions as in Fig. 3. 

2. Boj^O, Bi = 

The second case we consider is the case of a nuclear 
spin noise in a nondriven limit: £>o 7^ 0, B± = (u> = 
0). In this case the nuclear hyperfine field provides a 
flustuating field sampled by the electronic spin. Thus 
the total field will be given by a sum of the external 
field, Bo, and the AI Z term. In this case the Larmor line 
will acquire sidebands due to the nuclear hypcrfinc field. 
This is clearly seen from Eq. (8), which in the present 
case reduces to 

K+-(t) ~ |(1 + 8cos2AT)e- l2BoT , (10) 

since u>\ = |Aoj| = A giving Acj/oji = 1. Obviously, 
this expression provides a main peak at e = 2B$ and 
sidebands at s = 2B$ ± 2A : as is illustrated in Fig. 4 
(dotted). Similarly, the z-component of the spin-spin 
correlation function reduces to 

K z (t) ~ 3 - 2 cos 2Ar. (11) 

Thus, there are dips at e — ±2A, as seen in Fig. 4. Note 
that the amplitude of both and K z are independent 





6/60=0 








1/100 








V J 


k 


V 







e 



FIG. 4: Intensity of the normalised K + (e) + K z (s) with 
(solid) and without (dotted) rf field applied over the system. 

of the nuclear hyperfine field in this case, however the 
positions of the peaks obviously shift linearly with A. 

Finally we comment on the case with no external field, 
e.g. Bq,Bi = 0. Then, the impurity spin interact with 
the nuclear spin via the nuclear hyperfine field A. As- 
suming that the nuclear spin has a a very slow time de- 
pendent dynamics we can use the theory developed in 
the former section. Thus, under these circumstances the 
spin-spin correlation functions will reduce to 

K+-(t) - ^(1 + 8cos2Ar), (12a) 

K z {t) ~ 3-2cos2Ar. (12b) 

As expected, the z-component of the correlation function 
is unaffected by the absence of the field H , whereas the 
transverse component is translated to £0 — as Bq — > 
0. Hence, the spectrum for K z is the same as in Fig. 
4 (dotted), as well as the spectrum for , however, 
shifted to e = 0. 



III. CONCLUSION 

In conclusion we presented a theory for nuclear spin 
noise spectroscopy. We considered a specific case of an 
ESR STM where the nuclear spin dynamics is revealed in 
the tunnelling current noise. We argue that noise spec- 
troscopy is capable of detecting nuclear spin fluctuation 
via a hyperfine coupling to localized impurity electronic 
spin. We find that the spectrum of the noise is rich and 
depends sensitively on the rf frequency and detuning, Fig. 
3. The main features of the spectrum are i) the peak at 
the Larmor electronic frequency that acquires hyperfine 
split satellites. This part of the spectrum is coming from 
the transverse spin fluctuations. The ENDOR like phe- 
nomenon occurs in the noise spectrum where the nuclear 



spin flips due to the rf field with nuclear resonance fre- 
quency directly affect the electronic spin dynamics, ii) 
the peak at zero frequency (that realistically will always 
be obscured by the 1// noise) also acquires satellites at 
±A. This observation suggests that one measure the nu- 
clear spin dynamics even in the absence of external fields, 
as long as the hyperfine lines are outside the l/f noise 
peak. 



surface and a nuclear spin I (/ = 1) in the substrate. The 
spin Hamiltonian for this system can be written as 

H s = B S Z + AS I + B 1 I X cosujt + B 2 S Z coscut, (Al) 



with B ^> A^> Bi,B 2 , which yields the effective Hamil- 
tonian 
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APPENDIX A: DERIVATION OF THE SPIN 
SUSCEPTIBILITY 

In the paper we have been considering the interaction 
between an impurity spin S (S = 1/2) on the substrate 



J 



Hf = B S Z + ASJ Z + B x l x cos tot. (A2) 



1. K+- 



In order to see the effect of the electron spin flips 
on the susceptibility, we consider the correlation func- 
tion K+-(t 1 ,t 2 ) = tr {S+(ti)S(t 2 )p}, where S±(t 1 ) = 
cxv{iJ* 1 H<fdt}S ± eM-iCnfdt}, S+ = | T)eU I, 
S_ = I |) e (T |, and p = (| T) e (T I + I DeU I) ® U- Using 
the cyclic invariance under the trace gives 



K+-( tl ,t 2 ) = tr J (T le-^^'l T)U \e l ^ H ^ dt \ [) = e -«*o(t a - tl ) tr J (? | e -i/,?*f*| T)(i | e * |) 



_ e -»2S (t 2 -ti) tr 



(AI Z +B 1 I X cos ujt)dtrp e ift? (-AI Z +B 



\I X cosu)t)dt I 

y 



(A3) 



where TL c g = AS Z I Z + B\I X cosut. Transforming the system into the rotating reference frame gives 



ji e -i/ t * 1 2 (A7 z ±Si/ a! cosLut)dt _^ e -iuiI z T/2 e -i[AuI z /ivi±B 1 I x / (2ui)]u>i_t 



(A4) 



where the detuning parameter Aw = A — ui/2, oj\ = ^ Auj 2 + (Bi/2) 2 , and r = t 2 — t± have been introduced. We 
introduce some notation by defining Q± = (Aw7 2 ± BiI x /2)/lo\. In the algebra of the spin 1 operators we note that 
Il n = I z for positive integers n and I 2n+1 = I z for all non-negative integers n. The same relations hold for I x , i.e. 



j2n 



> 1, and I x n+1 — I x , n> 0. It is then easy to show that also Qj_ n = Q\, n > 1, and Q 



2n+l 



Q±, n>0. 



These rules for the algebra of the spin 1 operators give 



= 1 — A 2 (I — cos olt ) ± A sin CtT, 



(A5) 



where the operator A is either of I z , I x , or Q±, and a is a scalar. Using these identities along with the facts that 
tr n I z = tr n I x = tr n I z I x = tr n I z I x = tr n I 2 I z = 0, tr n I z I x = 1, tr n I% = tr n l 2 x = 2, and tr n l n = 3, we find that 
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the correlation function K + (t\,t2) = K + (r) is given by (recall that the imaginary part of the trace vanishes) 



K+-(r) = 




By tuning into resonance, i.e. Aui — which is given at u> — 2A and yields u>\ = Bi/2, it follows that 

K + -(t) = e - l2B « T \- + -cos2Ar + -cosSxt+ - cos(2 A + B^t + i cos(2.4 - B^r + cos - cos — 
[4 4 4 8 8 2 

+ i cos(2^ + B 1 /2)t + ]- cos(2A - B 1 /2)t - \ cos(A + B x )t - ]- cos(A - B{)t\. (A7) 



2. K z 



By means of the same approach we derive the z-component of the spin susceptibility K z (ti, £2) = tr {S z (ti)S z (t2)p}, 
where S z {h) = cxp{i f* 1 Uf dt}S z exp{-i f* 1 Ufdt} and 5 Z = (| TXT I - I l)U I)/ 2 - Using that [S z ,Hf] = we 
find, for the | |) (| |-component of p, 



AK z {t u t 2 ) = tr J (T [e^ 2 ^| TXT ^1 T) = tr J (T [e^ 2 TXT 1^ T) 



tr „ j 



rp e -i (AI Z +B 1 I X cos wt)dtji e i (AI Z +B 1 I X cos ut)dt 



J =tr„| 



■'} 



= x 1+5 (-) ) + 

cos(w + 2u)\)t 



cos ut + cos 2^it] 



+ 7 1 + 



1- 



Aw 



Auj 



1 



Aw 



Aw\2 



cos(w — 2wi)r 



+ 



V2o;i/ 



[cos — + cos wirj - 2 1 ( - — ) [cos r + cos r | 

2 V Wi / V 2wi / 2 



(j — 2^! 



1/Bi\4 w + 4wi w-4wi 

-I - — ) [cos(w + o>i)r + cos(w — loi)t + cos r + cos t\. 

2 V 2lj\ J 2 2 



(A8) 



At resonance Alo = 0we then have 



K z (t)alj=o = r flO + 6[coswt + cos2wit] + cos(w + 2wi)r + cos(w — 2u)\)t + 8[cos ^ + cosu>it] 
8 V 2 

r]) (A9) 



. r , , , . w + 4wi w - 4wi 
-4[cos(w + Wijr + cos(w — wi)t + cos r + cos t 
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